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Abstract. A method is introduced to derive resolution criteria for various a priori defined templates of brightness distribu¬ 
tion fitted to represent structures and objects in astronomical images. The method is used for deriving criteria for the mini¬ 
mum and maximum resolvable sizes of a template. The minimum resolvable size of a template is determined by the ratio of 
(SNR - 1 )/SNR, and the maximum detectable size is determined by the ratio of 1/5 NR. Application of these criteria is dis¬ 
cussed in connection to data from filled-aperture instruments and interferometers, accounting for different aperture shapes and 
the effects of Fourier sampling, tapering, apodization and visibility weighting. Specific resolution limits are derived for four 
different templates of brightness distribution: (1) two-dimensional Gaussian, (2) optically thin spherical shell, (3) disk of uni¬ 
form brightness, and (4) ring. The limiting resolution for these templates changes with SNR similarly to the quantum limit on 
resolution. Practical application of the resolution limits is discussed in two examples dealing with measurements of maximum 
brightness temperature in compact relativistic jets and assessments of morphology of young supernova remnants. 
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1. Introduction 

Accurate knowledge of resolution limits is one of the important 
aspects of astronomical observations and data analysis. There is 
a number of research fields in which a quantitative assessment 
of resolution limits is vital for interpretation of observations. 
Most celebrated examples include the star/galaxy separation in 
imaging surveys (e.g. Stoughton et al. 12002k measurements of 
stellar diameters (e.g. Tango & Davis l2002t . measurements of 
sizes and shapes of young supernova remnants (e.g. Marcaide 
et al. 119951 McDonald et al. 120011 Perez-Torres et al. 120021 
Bartel & Bietenholz 2003 1. and estimates of brightness temper¬ 
atures in extragalactic jets based on very long baseline interfer¬ 
ometry (VLBI) measurements of sizes of most compact regions 
(cores) of the jets (e.g. Moellenbrock et al. Il996i Lobanov et 
al. 120001 Horiuchi et al. l2004ll 

The Rayleigh (118791) resolution criterion and its several 
modifications are commonly considered as a measure of lim¬ 
iting resolution. The Rayleigh limit is expressed by the reso¬ 
lution distance R = 1.22 (/l/2ay), where A is the wavelength 
of light and ay = D/f is the ratio of the diffracting aper¬ 
ture, D, to the distance, /, between the aperture and the im¬ 
age plane. Strict application of the Rayleigh criterion is lim¬ 
ited to the specific case of resolving two point-like objects of 
comparable brightness. In terms of Fourier optics, resolution 
is better represented by the instrumental bandwidth that rep¬ 
resents the range of Fourier frequencies sampled by an imag¬ 
ing apparatus. Diffraction results in existence of the maximum 
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(cutoff) frequency q mdx that limits the resolution. Indeed, the 
Rayleigh limit has been reformulated, in Fourier optics, in 
terms of the Nyquist distance R = 7t/q mllx . Linear inversion 
methods applied to out-of-band extrapolation yield resolution 
limits smaller than the one prescribed by R. This provides the 
basis of the so-called “superresolution” techniques (Bertero & 
De Mol ll996l and references therein). 

In many practical cases, an observed brightness distri¬ 
bution has to be tested against a specific, extended pattern, 
for instance, a characteristic surface brightness distribution of 
a galaxy, stellar disk or supernova shell. A number of as¬ 
tronomical studies employ fitting image brightness distribu¬ 
tion by such pre-defined templates, with the definitions being 
based on physical models or previous observational evidence. 
Introducing a priori knowledge (or a postulate) of the bright¬ 
ness distribution of an emitting region is similar to using the 
out-of-band extrapolation (in which a definitve postulate is em¬ 
ployed about the shape of the emitting region). In both cases, 
the limiting resolution depends on the signal-to-noise (SNR) 
of detection. In contrast to the out-of-band extrapolation, the 
limiting resolution for a priori defined templates is determined 
ultimately by quantum fluctuations of light (Kolobov & Fabre 
120001 because the templates are represented by analytic func¬ 
tions that can be extended over an infinite range of frequencies 
in Fourier domain. 

An analytical criterion for calculating resolution limits can 
be derived for specific templates of brightness distribution 
(TBD, hereafter) from their respective Fourier visibility repre¬ 
sentations (see e.g. Thompson, Moran, Swensson 19861 Taylor, 
Carilli, Perlev ITW9t . Working in the Fourier plane has a spe- 
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cific advantage of separating positional information (contained 
in the visibility phase) from information about the extent of the 
emitting region (contained in visibility amplitudes). This has 
been recognized early on in radio interferometry and used for 
a variety of purposes, most notably for error estimates (see e.g. 
Fomalont ll999 t. 

Analysis of visibility distribution in Fourier domain is ap¬ 
plied in this paper to derive general resolution criteria and ob¬ 
taining resolution limits for specific TBD. The basic founda¬ 
tions of the method for deriving the resolution limits are sum¬ 
marized in Sect. [2] The general resolution criteria are intro¬ 
duced in Sect. 0 for determining the minimum and maximum 
resolvable size of a TBD. Several relevant TBD are described 
in Sect. 0 Their respective limiting resolutions are calculated 
in Sect.^and compared with the quantum limits on resolution. 
Applications of the resolution limits are outlined, in Sect. [6] 
with two examples dealing with studies of brightness tempera¬ 
ture in compact AGN and assessment of morphology of young 
supernovae. 


<7max is determined by the ratio D/A a b s of the diameter of the 
aperture to the observing wavelength. For interferometers, the 
diameter is substituted by the maximum baseline length, fi max , 
between individual elements. If q max is known, the Fourier 
transform of the brightness distribution in an astronomical im¬ 
age can be used for deriving rigorous resolution criteria for a 
TBD. To convert the result to angular units in the image plane, 
<?max has to be related to the size of the point spread function 
(PSF) or the synthesised beam, in the terminology adopted in 
radio interferometry. This can be achieved by considering the 
interferometric visibility sampling function (SF) — an effec¬ 
tive coverage of the Fourier plane, which is often also termed 
as the transfer function. The theoretical shape of the PSF is 
completely determined by the combination of the SF and the 
frequency range (q m j n , q m ax ). In practice, the shape and size of 
the PSF are further modified by aberrations, weighting and ta¬ 
pering (apodization) of data (e.g. Norton & Beer 119761 Briggs 
j1 9951 Briggs, Schwab & Sramek ll999 l. which also affects ac¬ 
tual resolution limits. These effects have to be taken into ac¬ 
count, thus reformulating q max in terms of the PSF. 


2. Basic definitions 

Consider an emitting region of finite size that has an integral 
flux density F and a brightness distribution 1(1, m). The bright¬ 
ness distribution 1(1, m) is measured in the image plane, which 
is described by a rectangular coordinate system (l, m). The cor¬ 
responding Fourier (or interferometric) visibility distribution 
V(u,v) is measured in Fourier plane (uv-plane), with the re¬ 
spective rectangular coordinate system (u, v). The two are re¬ 
lated via Fourier transform (T)'. 


V(u, v) = T 1(1, m) = 


/ OO /"»cx 

CO — Cx 


1(1, m) exp[-27T i(u l+v m)] dxdv. 


X OO /"*00 

I V(u,v)exp[2ni(ul+vm)]dxdy 

OO %J —CO 

The visibility distribution is complex-valued, with the visibility 
amplitude determining the limiting resolution and the visibility 
phase affecting ultimately the positional accuracy in the im¬ 
age plane. For a circularly symmetric brightness distribution, 
the relation between V(u, v) and 1(1, m) can be represented by 
Hankel transforms 


2.1. Sampling of the Fourier domain 

The shape of the PSF, B(l,m), depends on the SF (or uv- 
distribution) s(u, v) in the Fourier domain. The PSF is then, in 
the simplest case, B(l,m) = Ts(u,v)- The PSF is character¬ 
ized by its width at half power level (FWHM), which is often 
also termed the “half power beam width”(HPBW, hereafter). 
Generally, the two-dimensional shape of the PSF must be con¬ 
sidered, and resulting directional resolution limits should be 
derived. If this is not feasible, an equivalent circular HPBW, 
S c , can be used, which depends on the actual shape of the PSF 
or, more generally, on the shape of the SF. S c gives then a cor¬ 
responding size of a circular PSF that contains the same in¬ 
tegrated flux density as the elliptical (or generally speaking, 
non-circular PSF and SF). For an elliptical SF, 


V(q) 


f 


I(r)Jo(2nq r)rdr , 


J p»oo 

V(q)jQ(-2nqr)q d q, 

0 

where Jq is the zero order Bessel function and q — (u 1 + v 2 ) 1/2 
and r = (x 2 + y 2 )^ 2 are radial coordinates in the Fourier plane 
and image plane, respectively. For simplicity of derivation, cir¬ 
cular symmetry will be assumed throughout this paper, but the 
method of analysis is valid for the general case as well, which 
can be readily shown by expanding the variables r and q into 
their two dimensional representations (l, m) and (u, v). 

The relevant range (q m ; n , q maK ) of frequencies (coordinate 
dimensions) in the visibility plane is specific to each individual 
observational setup. In most of the practical cases, q m ; n —> 0 or 
<?min ^ <:/niax can be assumed. For filled-aperture instruments. 


S c =(ab) l/2 . 


(1) 


where a and b are the major and minor axes of the PSF. For a 
rectangular SF, 

S c = (t Tab) 1 ' 1 . (2) 

Equations m and □ can be applied to describe S c of most of 
the imaging instruments used in astronomy. 

The relation between the beam and the sampling function 
is even more complex if s(u, v) provides either incomplete or 
non-uniform coverage of the Fourier plane. The visibility sam¬ 
pling can then be characterized by the smallest and largest uv- 
distances (gmin, qwx) and by the density pa-(u,v) of the uv- 
samples (assuming g max » q mm )- The density is defined so that 
Pa-(u, v)d u dv gives the number of visibilities within the area 
given by (u + 0.5du , v + 0.5dv). The visibility sampling can be 
modified by a weighting function w(u, v) to emphasize certain 
fractions of the visibility distribution and optimize image sen¬ 
sitivity on selected spatial scales. With the weighting function 
applied. 


s(u, v) = Pcr(u, v)w(u, V') , 


(3) 
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giving the total number of visibility data points N(u,v ) = 
f f s(u , v)du dv. 

Depending on the form of s(u, v), an appropriate correc¬ 
tion factor S q must be applied to the equivalent (or directional) 
HPBW S c , which gives 

<7max = (SqScT 1 ■ (4) 

2.2. S q for a power law s{u, v) 

A power-law SF can be used to describe an ideal interferome¬ 
ter as well as a major fraction of single aperture instruments. 
The effect of a power-law SF on the HPBW can be calcu¬ 
lated analytically. Assuming circular symmetry yields s(q) = 
Po-iq)w(q). The density of samples is Po-(q) K q 2 - The func¬ 
tion w(q) is then given by c/ 1 , where /j is a power index decrib- 
ing the weighting. The choice of w(q) that yields s(q) = const is 
commonly known as “natural weighting”. Thus, P - 2, for the 
natural weighting, and the corresponding s(q) is given by 1 \(q), 
where I hr/j = H( 1 +q)H( 1 -q) and H is the unit Heaviside func¬ 
tion. The natural weighting increases sensitivity to extended 
emission at the expense of reducing the resolution of the im¬ 
age. The resolution is optimized by choosing w(q) = const 
(uniform weighting), for which the respective s(q) = Po-(q). 
This yields s(q) oc q 2 for the uniformly weighted synthesised 
interferometric aperture (similar to the result obtained for filled 
aperture telescopes). A general form of the sampling function 
is then 

s(q) = n(< 7 ) w(q)po-qdq = Il(q) cf~ 2 dq , (5) 

with p — 0 for the uniform weighting and p — 2 for the natural 

weighting. Thus, for a power law s(q) given by Eq. < 0 , 

*S q =2 1 ^ /2 , (6) 

and c/ max is related to the effective HPBW as follows: 

1 

<?max - 2 1 ~< 8 / 2 S c ’ ( ' 7 ' ) 

where a and b are the major and minor axes of the synthesised 

beam. For filled aperture instruments, p — 0 and A q = 1/2 
can be generally adopted. S q can be modified by apodization, 
depending on the shape and parameters of the apodizing func¬ 
tion. The range of p [0,2] corresponds to a range [-1,1] of the 
robustness parameter introduced by Briggs (TT9951 

In interferometers, the shape of the SF may deviate from 
the power-law due to several factors including sparsity of re¬ 
ceiving elements, non-optimal configurations, and reduced mu¬ 
tual visibility on very long baselines. A relative error cr q of *S q 
due to an incomplete SF of an interferometer can be approx¬ 
imated from the size of gaps Aq/q in the SF. For an interfer¬ 
ometer with N regularly spaced elements, cr q » (Aq/q) 2 = 
[(<7max/<7min) 3 ^ V - l] 2 (Lobanov 12003f ). The resulting cr q is a 
few per cent or less for most of the existing connected inter¬ 
ferometers. It is expected to be improved substantially in the 
future interferometric instruments such as ALMA and SKA. In 
present-day VLBI observations, cr q remains about 10%-15%. 

1 http://www.skatelescope.org/PDF/ska_memo38.pdf 


Thus, Eq. (J6) provides a relatively accurate assessment of N q . 
Numerical analysis of the specific instrumental SF should be 
employed whenever a better accuracy is required. 

2.3. Noise and SNR 

The Fourier transform has a fundamental property of preserv¬ 
ing the covariance of noise. It implies that the signal-to-noise 
ratio (SNR) measured in the image plane is the same as the 
one that would have been recovered if the same structure was 
fitted to the data in the visibility plane. Then SNR — F/cr n — 
V(0)/cr n , which is the ratio of the integral flux density to the 
variance of the additive white Gaussian noise with a power 
spectrum p(q) = cr 2 5(q). This property enables the derivation 
of resolution criteria from the Fourier transform of the image 
brightness distribution. It should be noted that the definition 
of SNR introduced above is different from another commonly 
used definition given by I(0)/<r n . 

The full SNR is applied for calculation of the resolution 
limits if the brightness distribution is fit by a single template 
or if the position of the template of interest is fixed, in a fit 
by multiple templates. If the brightness distribution is fit by 
multiple templates at variable positions, the SNR in the esti¬ 
mates of dn m and c/ res should be multiplied by a factor q = 

yl 1 - Cphase/^amp, where °"phase and ^amp is the noise in the vis¬ 
ibility phases and amplitudes, respectively. Typically, the noise 
is distributed evenly between the amplitudes and phases, and 
so cr phase — rr am p — cr n / a/ 2 and q — 1 / a/2. 

3. Resolution criteria 

Consider a circularly symmetric object in the image plane, 
which has an integral flux density F and a distribution of bright¬ 
ness I(r) that transforms into a visibility distribution V(q). The 
noise is given by cr n . The visibility distribution Viq) is normal¬ 
ized by F, implying V(O) = 1. The visibility distribution can be 
conveniently represented as 

V(q) = d'V s q, (8) 

where 'V s and d represent the shape and size of a TBD, respec¬ 
tively. 

3.1. Minimum resolvable size 

For a point-like brightness distribution, I(r) = 6(r), Viq) = 1 
at all q. For an object with extended brightness distribution, 
V(q) = 1 -&q), where £(q) depends on the size and shape of the 
object. If this deviation of the visibility distribution from unity 
can be detected in the data, then the brightness distribution I(r) 
is resolved. The function //(q) evaluated at q lmx can be related 
to the SNR of detection of the TBD. The resolution criterion is 
then given by £(q milx ) > <x n = 1 /SNR. This criterion implies 
that I(r) is resolved if 

V(<7max) < SN *Z 1 = US NR) . (9) 

Solving this equation for d , recovered from the expression for 
VTc/max), yields the limiting resolution dn m for a given template. 
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Following Eq. ©, V{q mdx ) = d\ n „ q max 'Vs, and the limiting res¬ 
olution is 


d I i r 


1 f m (S NR) 

<7 in ax *Vs 


( 10 ) 


The form 'Vs depends on the choice of the template. The ratio 


Klim = 


MS NR) 

"Vs 


( 11 ) 


provides a theoretical resolution factor independent from the 
size and shape of the PSF of the data. R\[ m depends only on the 
SNR and the shape of a TBD. 7? llm is a measure of the space- 
bandwidth product c = d q nrdx . Indeed 7?| lm = 2 cij m . This elu¬ 
cidates the connection between R lm and the Shannon number, 
S = d/R , which gives, for a spatial interval [-d/2, +d/ 2], the 
number of sampling points that are separated by the Nyquist 
distance R = n/q max (Toraldo di Francia fl 969t i. It follows im¬ 
mediately that 

Rl,m = nS\ im , (12) 

where the subscript in S implies that the Shannon number is 
evaluated at d\ lm . Recalling Eq. Q, the resolution limit can be 
rewritten as 

d fm =%, m S q S c . (13) 

The HPBW given by S c can represent both the FWHM of a 
non-circular PSF along a selected direction in the image plane 
or an equivalent of the FWHM of a circular PSF calculated 
using Eqs. {DO. The term *S q accounts for the effect of limited 
sampling of the Fourier domain (see Sect. 12.2> . 


and the largest resolvable size is then given by 
, _ SNR'^MS NR) o ^ 

^res — -&c 

y s 


(17) 


An equivalent of the SNR-dependent resolution factor 'R can 
now be defined by 


R 


res 


SNR l, °-Vf r (SNR) 


( 18 ) 


transforming Eq. o into the same form as Eq. O- By anal¬ 
ogy with Eq. ill 21 . 

Rrcs = n S res . (19) 


The limit d les signifies the largest size of a template that can 
be detected. At sizes larger than d res , sensitivity to extended 
emission becomes insufficient to detect the specific shape of 
brightness distribution described by the template. In case of fit¬ 
ting astronomical data by multiple templates (e.g. in “modelfit¬ 
ting” used in radio interferometry; see Pearson fl 9991 . d les gives 
an estimate of the size at which a template must be split into 
two templates, to maintain an adequate representation of the 
structure observed. The ratio d tes /du m defines an effective “res¬ 
olution dynamic range” (RDR), D, that expresses the range of 
template sizes that can be detected at a given SNR. The RDR 
depends only on the SNR of detection and the data weighting, 
since the terms related to HPBW, SF, and the shape of the tem¬ 
plate ('Vs) all cancel in the ratio, giving 

£> = . ( 20 ) 


3.2. Largest detectable size 

Visibility amplitude of a completely resolved template reaches 
cr n at a scale q ies such that q res < q mdx . This yields a resolution 
condition 

V(q res ) = cr n = ^ = f r (S NR). (14) 

The largest detectable size is obtained by using Eq. © to sep¬ 
arate different terms in V(q les ), which gives 


4. Relevant templates of brightness distribution 

There are several commonly used TBD for describing struc¬ 
tures observed in astronomical images (particularly in VLBI 
images). These include a two-dimensional Gaussian profile 
(“Gaussian component”), a disk of uniform brightness, an op¬ 
tically thin sphere (or shell of finite thickness), and a ring (see 
Pearson Gaussian components are commonly used to 

approximate and quantify structures recovered in VLBI images 
of relativistic jets and other objects. Shells, disks and rings are 
applied to analysis of images of supernovae and stellar objects. 


= J _mm ( i5) 

q Tes "V s 4.1. Gaussian component 


similar to Eq. {TOJ. The main difference from Eq. GD is that 
MS NR)/ < V S cannot now be used as an SNR-dependent resolu¬ 
tion factor separately from lcs , since the latter also depends on 
the SNR of the detection. Evidently, SNR = yjN{q m?lx ) / N(q les ) 
is required for a fully resolved template to be detected at q res . 
Here, N(q max ) and N(q ies ) are the numbers of visibility samples 
within < 7 max and q res , respectively. The value of q KS depends on 
the sampling function s(q), and it can be found by solving the 
equation N(q ) = J (| s(q)dq for the upper limit of integration. 

Alternatively, q res can be expressed as a combination of 
SNR and qwx, with the latter being conveniently related to the 
HPBW using Eq. Q. For instance, for a power-law s(q) given 
by Eq. ©. 


_ *7max 

qreS ~ SNR 1 ^ ’ 


(16) 


A circular Gaussian shape with the FWHM d is given by 


1(f) — 



-4 In 2 r 2 \ 

f j 


Co- 

T exp 




( 21 ) 


with Co- - 2 Vhi(2)//r. An elliptical Gaussian component 
I(r, (/>) can be described substituting d with the running FWHM 


d,/, = J[a cos((/> - i/0] 2 + [b sin(0 - i/r)] 2 , (22) 

where a, b are the major and minor axes of the component and 
(// is the position angle of the major axis. The visibility distri¬ 
bution corresponding to the Gaussian component is 


V(q) = exp 


-(nd# q) 1 'j 

4 In 2 j 


= exp 



(23) 
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An equivalent of circular FWHM can be obtained in this case 
by finding the corresponding diameter of a circular Gaussian 
which has the same area at the half-power level. This gives 
^fwhm = ' fab - 


4.2. Spherical shell of finite thickness 

A homogeneous shell with the outer radius r s and a thickness 
6 T r s = (1 - a)r s (with 0 < a < 1) is described by 


Hr) = 


fr s (r) - far,(r) 

0 


r < r s , 
r>r s 


(24) 



q 


with f z (r ) = Vz 2 - r 1 . Thus, a = 0 describes an infinitely thin 
shell and a —> 1 describes a filled sphere. The respective vis¬ 
ibility distribution is given by Hankel transform of Eq. Cl}, 
which yields 

V(q) =£ v 0 j , [sin(27r r s q) -2nr s q cos(27r r s q)- 

kn i q i (25) 

- sin(27r a r s q) + 2n a r s q cos(27r a r s q)\ , 

where £' v = 3/[r|(1 - a 3 )] is the normalization factor derived 
from the condition V(0) = 1. If the shell is described by its 
outer diameter d = 2 r s , then f z = yz 2 - d 2 /4. and 


Fig. 1. Visibility amplitudes V(q) of the templates plotted for 
the small-size (d\i m ) and large-size (d les ) limits. 

4.4. Ring 

In a ring of diameter d, the brightness is zero everywhere except 
on the circumference: 

Hr) = —5(r - d/2). (31) 

nd 

The respective visibility distribution is given by 

V(q) = Jo(ndq), (32) 


V ^4) = - ,, , ,, , [sin(zr dq) - ndqcos(ndq)- 

(1 -a i )n i d i q i (26) 

- sin(a ndq) + andq cos(a ndq)] , 


For small d (d < HPBW), V(q) can be approximated by the 
Taylor expansion 

v(q)=i -^[h^y d2qi+o(q4) - 


4.3. Disk of uniform brightness 

A uniformly bright disk of diameter d is described by 


(4/(nd 2 ), r < d, 
(0 r > d 


The corresponding visibility distribution is 


V(q) = 


2J\(ndq) 
nd q 


(29) 


where J\ is the Bessel function of the first kind. A useful ap¬ 
proximation of V(q) in the small-size limit is given by the 
Taylor expansion 


V(q) = 1 - 


(ndq ) 1 


+ o(q) 4 . 


(30) 


and it can be approximated conveniently using the Taylor ex¬ 
pansion of the Bessel function 


fotndq) - \+ o(q) A . (33) 

The visibility amplitudes of all five templates are compared 
in Fig.^l which illustrates the behavior of V(q) in the small-size 
and large-size limits (corresponding to the conditions used for 
deriving dn m and d rcs . respectively. It should be noted that the 
third order Taylor expansions used above approximate well the 
shape of V(q) only in the small-size limit. Thus, only d\i m de¬ 
rived from the Taylor approximations of V(q) are accurate. In 
the large-size limit, analytical expressions obtained with V(q) 
given by Eqs. and underestimate d les by up to 

40%. To remove this inaccuracy, numerical solutions must be 
obtained for c/ res of the shell, disk and ring templates. The nu¬ 
merical solutions can then be used to derive empirical correc¬ 
tion functions for the analytical expressions for d ies . These cor¬ 
rection functions are derived in AnncndixfAl The corrected d KS 
are accurate to within 0.01%. 

4.5. Noise term in Taylor expansions of V(q) 

The o(q 4 ) term in the Taylor expansions used for approximating 
V(q) must be expressed in terms of cr n or SNR, in order to be 
included into the derivation of d\\ m . The o(q 4 ) term reflects the 
error of the Taylor approximation of V(q). Since q represents 
a position in the two-dimensional Fourier plane, the error is 
proportional to the area cr 2 of the region described by o(q 4 ), 
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modified by the uncertainty 1/(1 - cr n ) in the border of that 
region. Therefore, o(q 4 ) can be expressed as follows: 

r 2 1 

(34) 


o(q 4 ) 


1 


SNR(SNR-l) 


As was discussed above, application of the Taylor expansions 
is restricted to the small-size limit, and the representation of 
o(q 4 ) given by Eq. can only be used for estimating d\\ m . 
Estimates of d les derived with this representation would diverge 
at small SNR. 


5. Resolution limits 


The resolution criteria given by Eqs. <n and (1 1 41 can be applied 
to derive specific resolution limits for the templates described 
above. Analytical solutions can be found for these equations, 
following the method outlined in Sect. 13.1H3.2I If analytical 
expressions for d\\ m and c/ res are computationally inconvenient 
or cannot be found, the limits can be obtained numerically. 
Numerical solutions should then be obtained for an effective 
HPBW = 1, thus requiring 


*7max 


1 1 

V ^ 2 1 -/ J /2 


(35) 


to be set for estimating <4m- For d KS 

_ j_i_L_ ™ 

qies ; 

should be set. Equations and <U3 can be then solved nu¬ 
merically for d\i m and d rcs expressed in units of HPBW. The 
corresponding resolution factors are obtained from 


'ftlim - d]im 22 -P /2 ’ ^ res _ ^ res 22-/S/2 ' ^7) 

These numerical solutions can be calculated for a relevant 
range of SNR, and the dependences of <4 rn and d ies on SNR 
can be established. 

This section and Appendix[A]summarize the analytical res¬ 
olution limits obtained for the templates described in Sect. 0] 
For each template, general forms of d\i m and d les and their re¬ 
spective 7?| lrn and 7? re s are listed. In all specific examples, a rect¬ 
angular, power-law SF given by Eq. 0 is adopted, with S c 
given by Eq. 0 and the SF correction factor N q described by 
Eq. 0. The uniform weighting (/3 - 0) is assumed in all ex¬ 
amples. This corresponds, with a high degree of accuracy, to 
astronomical images with a Gaussian PSF. 


5.1. Gaussian component 


d\im — 


n 


In 2 In 


/ SNR X 
\SNR-l). 


1/2 


Klin = - 
n 


In 2 In 


I SNR \ 
V S NR - 1 / 


S q S c , 

1/2 


For a rectangular, power-law SF, this corresponds to 


*4m — 


2 2-y3/2 


n a b In 2 In 


/ SNR \ 
\SNR- 1/J 


1/2 


(38) 

(39) 

(40) 


For the uniform weighting, Eq. TO yields a “benchmark” 
value of d hm = 1 HPBW at SNR = 4. 


5.1.1. Resolution limit in the image domain 


For the specific case of a Gaussian template fitted to an im¬ 
age with a Gaussian PSF, d\i m can be derived in the image 
plane as well. Consider the maximum difference between the 
Gaussian PSF, 4W> and a Gaussian template, /f(r), convolved 
with the PSF. In this case, the deconvolved size, <4, of the 
template can be conveniently expressed in units of FWHM of 
the PSF, thus setting the HPBW db = 1. This corresponds 

to a Gaussian with a FWHM d{ = 1 + dq fitted directly to 

the image. Normalizing the peak flux density and the PSF to 
/f(0) = 4(0) = 1 gives SNR = l/cr n . This normalization can 
be used for all df smaller than, or comparable with, the HPBW. 
The PSF is then described by 

4W = C £r exp(-41n2yj (41) 


and the template is given by 


4(0 = Co-exp 


-4 In 2 


1 +d\, 


(42) 


The template will be resolved in the image if 4-4 > cr n at 
least at one point. The maximum difference between 4 and 4 is 
realised at r — 1 / V2, and so the resolution condition becomes 


/ l \ !/(!+</?) 

*- 4 - i 



O- n _ 1 

C a ~ Co-SNR ' 


(43) 


Solving for d\i m = df gives the 


r/1 j in — 


ln(l/4) 


In (crJCr + 1/4) 


1/2 


ln(4) 


ln( 


4 CgSNR 
4+CV SNR 


) 


- 1 


1/2 

(44) 


measured in units of HPBW. This is the deconvolved size of a 
Gaussian template, and it can be compared with results from 
Eqs. H. For SNR = 4, Eq. I® yields d ]im = 1 HPBW, 
similarly to Eq. 6S3 applied to uniformly weighted interfero¬ 
metric data. The two limiting resolutions are essentially iden¬ 
tical at SNR > 4 (Fig 0. For SNR < 3, the asymptotes 
of Eqs. I 0D > and ®> are different, and Eq. ®> diverges at 
SNR = 4/3. This should be expected, remembering that the 
normalization used for 4 and 4 can only be used for df (and, 
consequently, <4; m ) smaller than the HPBW (which necessarily 
implies SNR > 4). 

For d\i m given by Eq. (1441 . the corresponding limiting size 
of a Gaussian template convolved with the PSF is given by 


d\u 


ln(4) 


ln( 4SNR ) 
' 4 +SNR ’ 


1/2 


(45) 


which is always larger than the HPBW. 


5.1.2. Maximum resolvable size and the RDR 

<4es = ^ [ln2 5/VR 2/(3 -^ln(5AR)] l/2 5 q S c , (46) 

K es = ^ [in 2 SNR 2/(3_/3) In (S /VR)]’ /2 . (47) 
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Fig. 2. Resolution limits derived for a Gaussian template in the 
Fourier plane (solid line) and image plane (dashed line). The 
resolution limit in the Fourier domain is derived for a rectan¬ 
gular SF and the uniform weighting. The resolution limit in the 
image domain is derived for a Gaussian PSF. At SNR > 4, both 
limits have the same dependence dn m oc SNR- 112 . At low SNR, 
the limit derived in the image plane diverges more rapidly be¬ 
cause the normalization 4(0) = 4(0) = 1 is valid only for 
d f < HPBW. 


The maximum resolvable size for a rectangular, power-law SF 
is given by 

92-jS/2 y /2 

d KS = —— [nab In 2 SNR 2 /®-® ln(SAR)] . (48) 


At very high SNR and for poor Fourier domain coverages, d les 
estimated from Eq. © can exceed the largest detectable scale 
cZids = l/<7min- In this case, di& should be used rather than c/ max . 
The respective RDR obtained from the ratio of d KS /du m given 
by Eqs. J46I and is then 


D = 


SNR 2 /®-® In (SNR)/ In ( SNR ) 
\SNR- 1/ 


1/2 


(49) 


The RDR is plotted in Fig.0for different values of /i. For D < 
1, no estimates of size can be made. The SNR levels at which 
D = 1 are 1.86, 1.82, and 1.70, for/? = 0,1,2, respectively. At 
SNR > 10, the RDR is oc SNR 016 , SNR 0M , and SNR 109 for 
P — 0,1,2, respectively. 


5.2. Quantum limits on resolution 

The dependence of <4 m on SNR can be compared with the ul¬ 
timate limit on resolution due to quantum fluctuations of light 
(Kolobov & Fabre 5000l . The effective quantum resolution dis¬ 
tance R q — RS / (Q+ 1), where Q + 1 is an effective number of 
degrees of freedom (NDF) in the template (Bertero & De Mol 
mm . The Rayleigh limit R = n/q mdx = n 3/2 2 X ~^ 2 yfab, and 
this yields the quantum limit of resolution 

R q = n N 2 2'-mJ_ (50) 

Q +1 

expressed in units of HPBW. The values of S and Q + 1 are 
calculated for each SNR and its corresponding i4 m . R q repre¬ 
sents then the absolute, quantum limit for resolving an emitting 


Fig. 3. Resolution dynamic range (D = d rcs /di mi ) derived for 
a Gaussian TBD, a rectangular SF and different values of [i. 
At high SNR and for poor Fourier domain coverages, d res may 
exceed <4i s and should be then substituted by the latter. Valid 
size estimates are only possible for D > 1. This corresponds to 
SNR > 1.86, 1.82, and 1.70, for /? = 0,1,2, respectively. 



Fig. 4. Resolution limit <4 m for a Gaussian template compared 
to the quantum limit R q of resolution and the superresolution 
limit for out-of-band extrapolation. At SNR > 3, the quan¬ 
tum limit and dn m are both oc SNR 1/2 and d\i m ~ 1.5 R q . The 
superresolution limit approaching oc SNR-° l at SNR > 10. 
SNR = 4.76 marks the lowest limit at which superresolution 
can be applied, and SNR = 20 gives the lowest limit for super¬ 
resolving a Gaussian brightness profile. 


region of finite extent and arbitrary shape detected at a given 
SNR. The same equation evaluated at S - n :,/2 2 ( A 2 t /2 gives 
the SNR dependence of the superresolution limit for out-of- 
band extrapolation. The NDF is evaluated following Kolobov 
& Fabre (120001 and Frieden Awn in the quantum limit and 
following Bertero & De Mol (1996)1 in the classical limit that 
applies for superresolution. 

The quantum and superresolution limits are compared to 
c/tim in Fig. 13 The quantum limit decreases at the same rate as 
d\i m , approaching asymptotically to an SNR- ]/2 proportionality 
at SNR > 3. This result agrees well with the theoretical limit 
on photon number fluctuation in an interferometer given by 
An « tr 1 / 2 = SNR- 1/2 (Forward 1 19781 . Hariharan & Sanders 
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fl996t . At SNR > 3, r/iim « 1.5 R q , thus connecting the reso¬ 
lution criterion for a Gaussian template to the quantum limit. 
Compared to the quantum limit and d\i m . the superresolution 
limit decreases much slower, approaching an S NR 0J depen¬ 
dence at SNR > 10. The difference between the two asymp¬ 
totes results from different strength of assumptions made about 
the extent and shape of the brightness distribution. The super¬ 
resolution is not feasible at SNR < 4.76 for an object of arbi¬ 
trary shape. For a Gaussian (or nearly Gaussian) brightness dis¬ 
tribution, application of the superresolution technique is only 
feasible at SNR >20. 

5.3. Comparison of resolution limits for different 
templates 

Resolution factors and limits for a spherical shell, uniform 
disk and ring templates described in Sect. 0] are derived in 
Appendix |A] Resolution factors 7?ii m and R res are plotted in 
Fig.0 and compared in Tables ^| an d 0 for different templates 
and SNR levels. The Gaussian template has the smallest ‘Rii m at 
all SNR. At SNR > 4, R hm oc SNR- 1 ' 2 for all templates. The 
picture is different for 7? res . At SNR —» oo, the asymptotes are 
R les x S NR oa for a Gaussian template and R ms oc S NR 033 for 
the other templates. This difference is caused by the presence 
of zeros in the V(q ) of the shell, disk and ring templates (see 
Fig.n). The Gaussian template reaches zero at infinity, and thus 
it can be better detected at large sizes and large SNR, compared 
to the other three templates. 


Table 1. R\[ m for different templates 


Rlim 

1 

3 

5 

SNR 

10 

10 2 

10 3 

Gaussian 

OO 

0.337 

0.250 

0.172 

0.053 

0.016 

Sphere 

oo 

0.712 

0.503 

0.335 

0.101 

0.032 

Shell 

oo 

0.551 

0.390 

0.260 

0.078 

0.024 

Disk 

oo 

0.637 

0.450 

0.300 

0.090 

0.029 

Ring 

oo 

0.450 

0.318 

0.212 

0.064 

0.020 


Table 2. R les for different templates 


R\im 

1 

3 

5 

SNR 

10 

10 2 

10 3 

Gaussian 

0 

0.801 

1.150 

1.733 

5.279 

13.930 

Sphere 

0 

1.396 

1.915 

2.701 

5.540 

14.171 

Shell 

0 

1.046 

1.413 

1.956 

4.595 

9.913 

Disk 

0 

1.231 

1.677 

2.345 

5.592 

12.182 

Ring 

0 

0.831 

1.111 

1.521 

3.524 

7.648 



Fig. 5. Resolution factors 7?| lm and 7? res for different templates. 
At a given SNR, the Gaussian template has the smallest 7?ii m . 
All templates have similar asymptotics 7?| lm oc SNR 1/2 at 
SNR > 4. For the largest detectable size and SNR —> oo, 
the asymptotes are 7? res oc S NR 0 4 for a Gaussian template and 
7?, cs oc SNR 0,33 for the other templates. The steeper rise of 7? res 
of the Gaussian template compared to the other templates re¬ 
sults from the presence of zeros in the V(q) of the shell, disk, 
and ring (see Fig.0. 


6. Applications 

6.1. Brightness temperature estimates from VLB I data 


Brightness temperature, 7k, of non-thermal continuum emis¬ 
sion is commonly used as an indicator of physical condi¬ 
tions in the emitting material. In relativistic jets, the inverse- 
Compton limit of 7 b ~ 10 12 K (Kellermann & Paulini-Toth 
1969 Kellermann 2002) is often reached and exceeded, im¬ 
plying angular dimensions smaller than the FWHM of the re¬ 
solving PSF. In this case, the resolution criteria described by 
7?i lm and d\ lm can be applied to analyze the size of a template 
fitted to interferometric visibilities or to image brightness dis¬ 
tribution. Comparison of c/| lm with the size obtained from the fit 
by a template can be used to separate valid measurements from 
upper limits. 

Elliptical or circular Gaussian components are applied rou¬ 
tinely for estimating sizes of emitting regions. The PSF is often 
also represented by an elliptical Gaussian shape. In the most 
general case an elliptical Gaussian component is described by 
its axes 0 ma j, (7 mln and the position angle of the major axis t //. 
Correspondingly, an elliptical PSF is given by (£> ma j, ^min. 01- 
In this case, the resolution limit should be calculated for f9 lnllJ 
and 0 m j n . The FWHM of the PSF is represented by its width by 
measured along the direction of the axis for which the limit is 
calculated. The corresponding resolution limit is then obtained 
from Eq. (ETil 


Aim,* = 2 2 - /3/ \ 


[ ln2 ln 

SNR \ 

[ 7T ln 

[SNR- 1/ 


1/2 


(51) 


Whenever either of the two axes is smaller than its respective 
0iim,i/,, the size obtained from the fit by a Gaussian component 
should be treated as upper limit on the size of the emitting re¬ 
gion. 
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Fig. 6. Zero order product logarithm function Wq(£) calculated 
for a range of shell width parameter a. 



Fig. 7. Comparison of the resolution limits for a Gaussian and 
an infinitely thin spherical shell. The shaded area marks angular 
sizes for which a fit by a thin shell cannot be distinguished from 
a fit by a Gaussian of the same size. 


6.2. Fine structure of young supernovae 

In VLBI observations of early stages of supernova expansion, it 
is important to decide at which moment the observed brightness 
distribution can be identified with a shell-like structure, thus 
distinguishing it from another shape (for instance, a Gaussian 
or a uniform disk). Let us assume that the brightness distribu¬ 
tion of emission from a young supernova is fit by a spherical 
shell of size d s and thickness a. The resolution criteria can then 
be used to establish whether (1) the fit is different from a fit by 
another template of the same size and (2) the fit can be distin¬ 
guished from a fit by another template of arbitrary size. These 
two cases form the necessary and sufficient conditions for veri¬ 
fying the shell-like shape of the expanding remnant. A method 
for distinguishing between a shell and a Gaussian templates is 
described below. The same procedure can be used to derive cri¬ 
teria for distinguishing between any other pair of TBD. 

To distinguish between a shell and a Gaussian component 
of the same size, d* s , the condition 


(52) 



Fig. 8. Comparison of the resolution limits for a Gaussian and 
a filled sphere. The shaded area marks angular sizes for which 
a fit by a filled sphere cannot be distinguished from a fit by a 
Gaussian of the same size. 


must be satisfied, where both the Gaussian (V„) and shell (V s ) 
visibility amplitudes are normalized to unity at q = 0, as given 
by Eqs. and O- This condition applies only to t/ g s com¬ 
parable or smaller than the HPBW. Setting HPBW = 1, the 
visibility difference in Eq. is maximized at q SiS = 1 / ^[\n. 
Equation can then be rewritten as 


ex P^T^d) + 
Ti 



1 

1 - cr 


(SNR - l) 2 + 1 
SNRiSNR - 1) ’ 


(53) 


where dj = Jd^ ?s /(ab) is the relative size of a template mea¬ 
sured in units of HPBW, and A\^ are given by: 


16 In 2 

di =- 

n 



n 



Denote 


fsNR - 


(SNR - l) 2 + 1 
SNRiSNR - 1) ' 


Solving Eq. (|53ji for <5(i and recovering the angular size b 
yields the minimum size at which it is possible to distinguish 
between the spherical and Gaussian shapes 


d\ t 


im,gs 


a b 
SNR 


[d2^SNR + 'llWo(f)] 


1/2 


(54) 


where 


'h 

T eXP 


-^ftNR 

'll 


and Wo is the zero order product logarithm function described 
in Appendix IbI Figure [6] shows the dependence of Wo(0 on 
SNR and the shell width a. The relevant range of SNR = 
[1, oo] determines the range of ( (-1/e, 0). Within this range, 
an approximation 


W 0 (£) + (1 + —)(£ + -) 1/e - 1 
e e 


can be used. The resulting limits on resolution are shown in 
Figs-EHElfor a thin shell and filled sphere, respectively. In each 


/(Vs, V g ) = | max(V s - V g )|«“ = cr n = 1/SNR 
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Fig. 9. Comparison of the resolution limits for a Gaussian and 
an infinitely thin spherical shell. The shaded area marks angular 
sizes for which a fit by a thin shell cannot be distinguished from 
a fit by a Gaussian of arbitrary size. 



Fig. 10. Comparison of the resolution limits for a Gaussian and 
a filled sphere. The shaded area marks angular sizes for which 
a fit by a filled sphere cannot be distinguished from a fit by a 
Gaussian of arbitrary size. 


Appendix A: Resolution limits for selected 
templates 

This appendix summarizes analytical expressions derived for 
the resolution factors 7?| lln and 7? ICS and limiting sizes r/| lln and 
d rcs . The analytical expressions for K res and d Tes are accurate 
within 40%. Empirical correction factors k(SNR) are intro¬ 
duced to correct this inaccuracy. The correction reduces the 
errors to within 0.01% for SNR < 1000. At larger SNR, the 
correction factors become essentially constant. 


A. 1. Spherical shell of finite thickness 


d]\m — 


1 - a 3 


Klim = - 
n 


1 - a 5 2 (SNR - 1) 
~ 1 - a 3 5 


1/2 


s q s c . 

(A.l) 

1/2 

(A.2) 


1 - a 5 2(SNR — 1) 

For a rectangular, power-law SF and uniform weighting, the 
corresponding limiting resolution for a spherical shell becomes 



40 /1 — a 3 \ cib 


- or 5 I SNR - 1 ' 


(A.3) 


The two limiting cases of the spherical shell are described by 
the filled sphere (a = 0) and infinitely thin shell (a —> 1), 
which gives 


r^lim,sphere 


^lim,shell 


/ 40 

ab 

1 71 

SNR - 1 

( 24 

a b 


n SNR - 1 
The largest detectable size is given by 


d KS = - c HSNR) 
n 


5 1 


2 1 


1/2 


s q s c 


(A.4) 

(A.5) 


(A.6) 


2 [5 1 -or 3 ' 

Kes=-<f KSNR) -- - J , 

n [ 2 1 - a 3 

where 

®(SNR) = 5AR 1/(3 ~ ffl ( 5Af/; ~ M 
\ SNR ) 


(A.7) 


figure, the shaded area indicate sizes for which it is not pos¬ 
sible to distinguish between a Gaussian and a spherical shape 
of the same size. At large SNR levels this limit approaches the 
resolution limit for a spherical shape. 

For a Gaussian and a shell of arbitrary sizes d g , d s , the vis¬ 
ibility difference reaches a maximum at q gs < q mllx , which is 
determined by solving the equation df 2 (V g ,Vs)/dq = 0. For 
each given d s (or d g , reciprocally), the minimum SNR at which 
one template can be distinguished from another is given by 
2//(V s , V g ) evaluated at q gs . Below this SNR, it is not possi¬ 
ble to claim that the template fit is unique. Application of this 
method is exemplified in Figs. IhlllOlfor a thin shell and filled 
sphere matched against a Gaussian. Similar approach can be 
applied to any other pair of templates. 


A.2. Disk of uniform brightness 

The resolution limit is given by 

d lim — 


2 

Klim = - 
n 


SNR- 1 
2 


1/2 


5qS c , 

1/2 


SNR- 1 

The largest detectable size is determined by 
d KS = - yf2®(SNR)S q S c , 


K res = - V2 0 (SNR). 


(A.8) 

(A.9) 


(A. 10) 
(A.11) 






























A.3. Ring 

The resolution limit is given by 


Lobanov: Resolution limits in astronomical images 


11 


, _ 2 

^lim 

71 


l 


SNR- 1 


1/2 


s q s c 


Klim = " 
n 


1 


1/2 


SNR- 1 

The largest detectable size is described by 


d les = —0(SNR)SnS c , 

n 


Kes = -0(SNR). 

7T 


(A. 12) 


(A.13) 


(A. 14) 


(A.15) 



Fig.B.l. Zero order product logarithm function Wo plotted for 
a range of arguments that cover the small and the large argu¬ 
ment asymptotics. 


A.4. Empirical corrections forties and 

Analytical expressions for d les and 7? res must be corrected for 
the error of the Taylor expansion in the large-size limit. Thus, 
for the shell, disk, and ring templates 

^res,corr — K (SNR)<R ies , (A. 16) 

with the empirical correction factors k(SNR) given by 

k(SNR) = 1 lg (a 2 lg (SNR) + 1) + 

1 \ (A. 17) 

a 4 [lg(5/V/?)] 2 + l/ ' 

The coefficients a\-ci 4 are listed in Table IA. fl At SNR > 300, 
k(SNR) ~ const, for all templates. The limiting values of 
k(SNR) are given in Table lA.il 

Table A.l. Coefficients for k(SNR) 



Sphere 

Shell 

Disk 

Ring 

d\ 

8.95923 

6.30926 

6.87495 

3.24248 

d 2 

0.11622 

0.13358 

0.13955 

0.20988 

a 3 

2.23652 

1.34282 

1.69345 

0.82607 

d4 

0.12907 

0.22302 

0.17399 

0.31843 

k(S NR > 300) 

1.42 

1.28 

1.35 

1.20 


Appendix B: Product logarithm function 

The zero order product logarithm function Wo (Fig. ED pro¬ 
vides the principal solution for w in the equation wexp(w) = 
z{w). Wo is the principal value of the Lambert W-function 
(Corless et al. Q3I13, and it is essentially an extension of the 
logarithm function. Wo is real-valued for z > -1/e. 
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